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Abstract. A 0.242 Mq object that finally becomes a helium white dwarf is 
evolved from Roche lobe detachment down to very low luminosities. In doing 
so, we employ our stellar code to which we have added a set of routines that 
compute the effects due to gravitational settling, and chemical and thermal 
diffusion. Initial model is constructed by abstracting mass to a 1 AIq red giant 
branch model up to the moment at which the model begins to evolve bluewards. 
We find that element diffusion introduces noticeable changes in the internal 
structure of the star. In particular, models undergo three thermonuclear flashes 
instead of one flash as we found with the standard treatment. This fact has a 
large impact on the total mass fraction of hydrogen left in the star at entering 
the final cooling track. As a result, at late stages of evolution models with 
diffusion are characterized by a much smaller nuclear energy release, and they 
evolve significantly faster compared to those found with the standard treatment. 

We find that models in which diffusion is considered predict evolutionary 
ages for the white dwarf companion to the millisecond pulsar PSR B1855+09 in 
good agreement with the spin-down age of the pulsar. 



1. INTRODUCTION 

Low mass, helium white dwarfs (WDs) have drawn the attention of many re- 
searchers who have devoted a large effort to their study. Currently, it is well 
established that helium WDs should have been formed in a binary system be- 
cause an isolated star would reach such a configuration in a time-scale much 
longer than the present age of our Universe. Examples of recent studies on he- 
lium WD evolution are Althaus & Benvenuto (1997; 2000), Benvenuto & Althaus 
(1998), Hansen & Phinney (1998), Driebe et al. (1998) and Sarna, Antipova & 
Ergma (1999). 

A crucial issue in the study of helium WDs is the amount of hydrogen left 
atop the WD after the end of mass transfer episodes. Sarna et al. (1999) have 
presented detailed calculations of the pre - WD evolution during the mass trans- 
fer stage in a close binary system. They found that, after detachment, helium 
WDs are surrounded by massive hydrogen envelopes of 0.01 - 0.06 Mq with a 
surface hydrogen abundance by mass of Xjr = 0.35 - 0.5. Massive hydrogen 
envelopes have also been found by Driebe et al. (1998), who have mimicked the 
binary evolution by abstracting mass from a IMq model at an adequate rate. 
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Binary systems composed by a helium WD and a millisecond pulsar of- 
fer us a good opportunity to test the predictions of evolutionary calculations. 
Millisecond pulsars are thought to be recycled during the mass transfer stage 
(Bhattacharya & van den Heuvel 1991). When mass transfer ends, the pulsar 
begins to spin down. If we set zero ages at the end of such mass transfer, then 
the ages of the WD and the pulsar component should be the same. If the pulsar 
radiates as a dipole and its current period (P) is much longer than the initial 
one, then the spin-down of the pulsar is given by tpsn — P/2P, where P is the 
rate of change of P. 

Recently, van Kerkwijk et al. (2000) have optically detected the low-mass 
component of the PSR B1855+09 system and determined its effective tempera- 
ture to be 4800 ± 800 K. In addition, the mass of the WD is accurately known 
thanks to the measure of the Shapiro delay of the pulsar signal and results to 
be O.258^o oi6-^0- The a S e derived for the WD component results to be twice 
as long as the pulsar spin-down age (« 5Gyr) when using the models without 
diffusion of Driebe et al. (1998). 

Helium WD evolution in presence of element diffusion has been considered 
by Iben & Tutukov (1986), who computed the evolution of a 0.296Mq object 
throughout the entire evolutionary stages and found two thermonuclear flashes. 
Remarkably, they found that after flashes the amount of hydrogen in the envelope 
of the model was unable to support an appreciable amount of hydrogen burning 
reactions, enforcing a fast evolution for the resulting WD star. 

It is the purpose of this work to study the impact of element diffusion on 
the structure and evolution of a 0.242Mq helium WD star. This mass value is 
representative of the companion to PSR B1855+09. In doing so, two sequences 
are computed: one with and the other without diffusion. Thus, the resulting 
differences should be due to the allowance for such a process. 

2. EVOLUTIONARY RESULTS 

For the present work, we shall employ the same stellar code we used in our 
previous works on WD evolution (Althaus & Benvenuto 1997; Bcnvcnuto & 
Althaus 1998). For the computations presented in this paper, we have included 
gravitational settling, chemical and thermal diffusion, as described by Althaus 
& Benvenuto (2000), and we refer the reader to that paper for further details. 

In order to have a realistic starting model, we have abstracted mass to 
a IMq giant star up to the moment at which it begins to evolve bluewards. 
In this way we have simulated the mass transfer episode during Roche lobe 
overflow. From then on, the evolution of the model has been followed assuming a 
constant value for the stellar mass down to log L/Lq = — 5. The results we have 
calculated without diffusion arc in very good agreement with those presented by 
Driebe et al. (1998), thus in comparing the results we have found considering 
and neglecting diffusion we are indeed also comparing with the results of Driebe 
et al. 

The evolutionary tracks for the models with and without diffusion are shown 
in Fig. 1. While the standard treatment predicts the occurrence of one loop, 
models with diffusion undergo three loops. Each of these loops are due to unsta- 
ble nuclear burning at the bottom of the hydrogen-rich envelope, i.e., thermonu- 
clear flashes. These flashes are critical events that, to a large extent, determine 
the subsequent evolution of the object, even during the final WD cooling phase. 
The occurrence of three flashes in models with diffusion is in sharp contrast with 
the predictions of the standard treatment, and clearly indicate that element dif- 
fusion is by no means a negligible phenomenon in the context of such stars. 
After flashes there are somewhat long time periods for which both nuclear and 
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Figure 1. Evolutionary tracks for a 0.242 Mq object which finally reaches 
a helium WD structure considering and neglecting element diffusion. Solid 
line corresponds to models with diffusion, whereas short dashed line depicts 
standard models without diffusion. Notice that models with diffusion undergo 
three loop excursions to the red giant part of the HR diagram, whereas models 
without diffusion experience only one of such episodes. 

photon luminosities are high. These periods correspond to evolutionary stages 
during which the total amount of hydrogen in the star (M#/M*) is noticeably 
reduced, as it is apparent from Fig. 2. 

An appreciable hydrogen envelope reduction is also found after the only 
flash in the sequence of models without diffusion but the key difference between 
the two sets of models is related with the occurrence of nuclear burning during 
the final WD cooling phase. In models with diffusion and after flashes, the total 
mass fraction of hydrogen left is Mh/M* w 1.8 x 10~ 3 . However, in the case of 
models without diffusion, when the star enters the cooling track, the amount of 
hydrogen that remains is about three times as large as in the case with diffusion. 
As a direct consequence of such a large difference, models without diffusion burn 
~ 2/3 of its hydrogen content during the cooling track. In models with diffusion, 
nuclear reactions also occur, but at a much lower total rate. This fact has, as 
we shall show below, a large impact on the final age of the star. 
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Figure 2. The total mass fraction of hydrogen as a function of time for the 
0.242 Mq helium WD model considering and neglecting element diffusion. 
Solid line denotes models with diffusion and dotted line corresponds to models 
without diffusion. Ages corresponding to the occurrence of thermonuclear 
flashes are indicated by arrows. A noticeable hydrogen depletion is found 
immediately after flashes. Notice that models without diffusion enter the final 
cooling branch (t > lGyr) with a much higher hydrogen mass, which is burnt 
during this epoch. In models with diffusion, nuclear reactions also occur, but 
at a much lower total rate. This fact has a large impact on the final age of the 
star. 

From the results presented in Fig. 3 (upper panel) it is clear that models 
with diffusion undergo a much faster cooling compared with models without 
diffusion. In models without diffusion, nuclear reactions release an amount of 
energy very close to the radiated one. This strongly delays cooling to very long 
ages. On the contrary, diffusion leads to hydrogen envelopes thin enough so as to 
prevent nuclear reactions from dominating the energetics of the star. Thus, the 
star has to extract energy from the relic thermal content, forcing a much faster 
cooling (in agreement with results presented in Benvenuto & Althaus 1998). 
To put this assertion on a more quantitative basis, we note that in reaching 
log L/Lq = —4 models with diffusion require 5.3 Gry, whilst models without 
diffusion need 13 Gyr. In view of the present results, the controversy found by 
van Kerkwijk et al. (2000) about the large discrepancy between the ages of the 
pulsar PSR B1855+09 and its WD companion is strongly alleviated. 

3. Discussion and Conclusions 

From the results presented in the previous section, it is clear that element dif- 
fusion introduces important changes in the internal structure of the star. Our 
present detailed results indicate that, contrary to Driebe et al. (1998) and 
Schonberner et al. (2000) assertions, nuclear reactions are not an important 
ingredient in determining the age of the star during the final cooling evolution. 
This conclusion is valid at least for the stellar mass value analyzed in this paper. 
Although the hydrogen mass fraction at the end of the computed stages is very 
similar in both sequences with and without diffusion, there is a large difference in 
the total time spent in reaching such a low luminosity in each case. While stan- 
dard models without diffusion spend about 30 Gyr in reaching logL/LQ = — 5, 
models with diffusion get such a luminosity in about 15 Gyr. 
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Figure 3. Photon and nuclear luminosity vs. age relationship for models 
considering (solid lines) and neglecting diffusion (dashed lines). Clearly, mod- 
els with diffusion evolve much more rapidly because the hydrogen envelope 
left is so thin that it is unable to support a appreciable amount of nuclear 
burning (lower panel). This is in contrast with standard results without dif- 
fusion for which hydrogen burning is the main energy source even at very low 
luminosities. 

As stated in the introductory section, previous models available in the liter- 
ature (those in which hydrogen burning is dominant during most of the cooling 
stages) lead to a strong discrepancy between the ages of PSR B1855+09 and 
its helium WD companion. These ages should be, on very general grounds, the 
same. The results presented in this paper strongly indicate that such a discrep- 
ancy is a direct consequence of ignoring element diffusion in the stellar model 
calculations. If we include it properly, ages naturally come into a nice agreement. 
Consequently, we do not have to invoke any ad-hoc mass loss as proposed by 
Schonberner et al. (2000) or exotic mechanisms to account for the observational 
situation of the binary system containing PSR B1855+09. 

It is worth noting that the structure of the models belonging to the final 
cooling branch we have found in the present paper for the 0.242 Mq helium WD 
strongly resemble those we have assumed in our previous work on helium WDs 
with hydrogen envelopes (Benvenuto & Althaus 1998). In fact, if we assume 
Mh/M* = 2 x 10~ 3 the agreement with the present results is excellent, not only 
regarding the cooling evolution but also with respect to the nuclear activity. In 
our opinion, the present detailed results largely justify the assumptions we have 
made in our quoted previous paper. 

We want to emphasize that our results strongly resemble those presented 
sometime ago by Iben & Tutukov (1986) (IT86) for the case of a O.296M 
model. They found two flashes and suspected that the second one was due 
to diffusion, but because they did not compute models without diffusion, they 
were not conclusive about this point. Here we confirm that diffusion is the 
reason for the occurrence of additional flashes. IT86 also found that the final 
amount of hydrogen left in the WD is small enough to ultimately enforce a 
fast cooling. Although their model is more massive than ours, which prevents 
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from a quantitative comparison between both set of calculations, the role of 
hydrogen burning at the final cooling stages is the same we have found. There 
is nevertheless an important difference between the work of IT86 and our study. 
In fact, IT86 considered that shortly after the second flash, when the model 
gets its maximum radius, it again overfills its Roche lobe and losses a tiny part 
of its stellar mass. Obviously, this mass belongs to the hydrogen envelope and 
such a loss has a large impact on the final value of Mfj/M» in the WD. This is 
ultimately responsible for the fast cooling occurring at late stages of evolution. 
The point to emphasize here is that we did not consider such a mass transfer 
during our calculations, and some words about this point are in order. If such 
a transfer actually occurs, Mh/M* should drop below the values we presented 
above (eventually to zero). This would have an impact on the global structure 
of the WD. Such a situation may be considered in the frame of the models 
presented in Benvenuto & Althaus (1998). Apart from the structural differences 
we do not expect an important change in the time-scale of cooling. Notice that 
in the models presented in Benvenuto & Althaus, the time spent in reaching a 
given luminosity value is rather independent from the exact value of Mh /M* (if 
it is in the range of I x 10~ 3 to zero). Thus, the occurence of such a late mass 
transfer should not modify the conclusions about the age of the final WD model 
and its better agreement with the observations related to the ages of the PSR 
B1855+09 system. 
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